Multiple pre-edge structures in Cu K-edge of high T c cuprates 
revealed by high resolution x-ray absorption spectroscopy 
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Using high resolution x-ray absorption spectroscopy and state-of-the-art electronic structure cal- 
culations we demonstrate that the pre-edge region at the Cu K-edge of high T c cuprates is composed 
of several excitations invisible in standard X-ray absorption spectra. We consider in detail the case 
of Ca2-a:Cu02Cl2 and show that the many pre-edge excitations (two for c-axis polarization, four 
for in-plane polarization and out-of-plane incident X-ray momentum) are dominated by off-site 
transitions and intersite hybridization. This demonstrates the relevance of approaches beyond the 
single-site model for the description of the pre-edges of correlated materials. Finally, we show the 
occurrence of a doubling of the main edge peak that is most visible when the polarization is along 
the c-axis. This doubling, that has not been seen in any previous absorption data in cuprates, 
is not reproduced by first principles calculations. We suggest that this peak is due to many-body 
charge-transfer excitations, while all the other visible far-edge structures are single particle in origin. 
Our work indicates that previous interpretations of the Cu K-edge X-ray absorption spectra in high 
T c cuprates can be profitably reconsidered. 

PACS numbers: 74.70.Ad, 74.25.Kc, 74.25. Jb, 71.15.Mb 



I. INTRODUCTION 



The interpretation of X-ray absorption spectroscopy 
(XAS) at the Cu K-edge in high-T c cuprates, and es- 
pecially of the pre-edge features, remains a challenging 
issue due to the strongly correlated nature of these mate- 
rials. Understanding the nature of the pre-edge features 
is relevant since these excitations are low in energy and 
probe the local environment of the absorbing atom. In 
several recent works [ll-Ql the presence of off-site excita- 
tions in transition metal compounds like NiO, cuprates 
and cobaltates was demonstrated. Off-site excitations 
can occur by a direct transition of a Is core electron to 
empty electronic states of the nearest neighboring atoms 
0. This mechanism has a low intensity (~ 1% of the edge 
jump) since the p-states of the absorbing atoms are not 
involved. More frequently the off-site excitations are me- 
diated by the on-site p-states of the absorbing atom that 
are very extended and therefore hybridize with neighbor- 
ing sites. When excitations due to intersite Cu (absorb- 
ing atom) 4p - Cu (other atom) 3d hybridization occur, 
the information is very relevant since these features are 
typically weakly affected by core-hole attraction and thus 
they probe the excitation in the absence of a core-hole in 
the final state. In some cases off-site transitions can be 
used to measure the charge transfer gap, as was recently 



suggested for NiO [|| . Therefore, a full understanding of 
pre-edge features and of their possible non-local nature 
leads to the knowledge of the hybridization mechanism 
between different orbitals in the low energy region (up to 
some eV from the Fermi level) and of the position of the 
upper Hubbard band in the absence of a core-hole in the 
final state [2( . 

In the case of high T c cuprates, and of correlated ma- 
terials in general, the occurrence of many-body charge 
transfer excitations in the edge region has been widely 
debated 0-0] • Recently, using first principles calcula- 
tions, we suggested Q that edge and far-edge structures 
previously attributed to charge transfer effects were ac- 
tually single particle in origin. Thus most of the multi- 
determinant effects seen in photoemission seem to be 
suppressed in absorption, probably by the presence of 
a final-state in the XAS cross-section. 

In this work we present high resolution angular re- 
solved partial-fluorescence yield (PFY) x-ray absorption 
measurements and first principle electronic structure cal- 
culations of the recently synthesized Ca2-a;Cu02Cl2 oxy- 
chlorides [TtJ ■ These materials form a new family of high 
Tc superconductors with a critical temperature peaking 
at 43 K at optimal nominal doping level (after post- 
annealing). The oxychlorides are appropriate for in- 
vestigating the electronic properties of the Cu02 plane. 
They are, for example, well suited for calculations be- 



2 




FIG. 1: (color online) Left: Crystal structure of 
Ca2-a;Cu02Cl2 0. Right: 14 atoms monoclinic supercell 
used in the calculation. 



cause of lesser structural disorder and the absence of rare 
earth atoms. They also present peculiar features com- 
pared to the cuprates such as a checker-board pattern 
of charge inhomogeneities at odds with the stripes 
in the cuprates. High resolution XAS allows us to iden- 
tify new pre-edge electronic excitations and a doubling 
of the main edge that are not visible in previous XAS 
measurements in high T c cuprates. Importantly we find 
that most of these excitations are single particle in na- 
ture. After identifying the single particle excitations and 
addressing the role of many-body charge-transfer effects 
in the XAS spectra, we come to the conclusion that the 
interpretation of Cu K-edge XAS in cuprates should be 
revisited. 

Section [TT] presents the experimental and theoretical 
details, and draws a comparison between the measured 
PFY-XAS spectra of Cai. 8 Cu0 2 Cl 2 and of La 2 Cu0 4 . In 
section IIII1 we confront the experimental data with the 
calculated cross section and identify the single particle 
peaks using the projected density of states on selected 
electronic states in the presence of a core-hole in the Is 
state. Finally we discuss the relevance of our results for 
the understanding of the empty state electronic structure 
in high T c cuprates. 



II. TECHNICAL DETAILS 
A. Experiment 

The experiment was carried out at SPring-8 on the 
inelastic x-ray scattering (IXS) beamlinc BL12XU. The 
incident energy was selected by a Si(lll) double crys- 
tal coupled to a high-resolution four-bounce Si(400) 
monochromator providing an estimated bandwidth of 220 
meV. We have used well characterized superconducting 
single crystals of Ca 2 -a;Cu02Cl2 (x — 0.2)[l7|. To pre- 
vent degradation of the samples known to be highly hy- 



groscopic, the crystals were covered with a grease suitable 
for low temperature measurements. The superconduct- 
ing temperature of 35 K measured by SQUID confirms 
that the intrinsic quality was preserved. The sample ori- 
entation was verified by x-ray diffraction prior to mea- 
surements. 

The absorption spectra were recorded in the partial 
fluorescence yield (PFY) mode using the beamline IXS 
spectrometer. The PFY-XAS consists of measuring the 
intensity variations of the Cu Kq (2p — > Is) emission 
line while scanning the incident energy through the Cu 
K-edge. With respect to conventional XAS, the PFY 
method yields an improved intrinsic resolution as the 
lifetime broadening effect in the PFY final state (with 
a 2p core hole) is considerably smaller than that of stan- 
dard XAS at the K-edge (with a Is core hole). The 
spectrometer was equipped with a 2 m bending radius 
Si(444) analyzer working at Bragg angle of 79.3° at the 
Cu Ka energy (8048 eV) and a high flux silicon drift de- 
tector. The total energy resolution was estimated at 250 
meV from the width of the elastic line. The sample was 
mounted on a goniometer head with the c axis (normal 
to the surface) in the horizontal plane. The PFY ab- 
sorption spectra were acquired with the polarization of 
the incident photons set either parallel or perpendicular 
to the c axis upon rotating the sample around a vertical 
axis. The spectra were corrected for self-absorption effect 
using the FLUO code developed by D. Haskel (Advanced 
Photon Source, Argonne). 



B. Calculation 

Ab initio calculations on Ca2-a;Cu02Cl2 are performed 
using spin polarized density functional theory (DFT) 
with PBE exchange-correlation [l2| . We cross-check our 
results with the DFT+U approximation, with U = 9.6 
eV as used in other cuprates @. The calculation are 
performed with the XSPECTRA code @, [lj] available 
under the GPL licence in the quantum-espresso package 
[l4j . We use Ultrasoft pseudopotentials [l5[ for all atomic 
species. The Is core-hole is included in the pseudopoten- 
tial of the absorbing atom. For the Cu pseudopotential 
we use non-linear core corrections [16JI and do not include 
semicore states. It has been shown [8| that the inclusion 
of the copper semicore states in the valence states is not 
necessary for an accurate description of XAS in cuprates. 

For the XAS calculation we use a monoclinic supercell 
containing one CuO plane and 14 atoms with experimen- 
tal lattice parameters [l7[ (see Fig. [T] ) . 

The magnetic state of Ca2Cu02Cl2 has not been yet 
resolved, but for simplicity we assume a collinear antifcr- 
romagnetism, similar to other cuprates. 

The wavefunctions are expanded on a plane wave basis 
with a 30 Ry cutoff for the wavefunctions, and 400 Ry for 
the charge density. A uniform Monkhorst-Pack 6 x 6 x 2 k- 
points grid (with respect to the antiferromagnetric unit 
cell) is used for the charge density calculation, with a 
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FIG. 2: Comparison between the PFY-XAS Cu K-edge spec- 
tra of Ca2-a;Cu02Cl2 (solid lines) and La2CuC>4 (dashed 
lines) for polarizations e || c and e || ab respectively in the 
top and bottom panels. In the case of Ca2-iCu02Cl2 the 
scattering geometries considered are: (i) k parallel to the ab 
plane and e parallel to the c-axis, labeled shortly e || c and (ii) 
k parallel to the c-axis and e parallel to the ab plane, labeled 
e || ab. In the case of La2Cu04 the geometries are slightly 
different, as reported in ref. p[. 

Methfessel-Paxton smearing of 0.05 Ry. The hole doping 
is treated by adding a compensating charge background, 
leading at x — 0.2 to a nonmagnetic state in agreement 
with experiments. 

The XAS cross section is calculated in the dipole ap- 
proximation using an ultrasoft-based continued fraction 
approach Q and PAW reconstruction [18]. Two paw pro- 
jectors are used in the absorbing atom pseudopotential. 
A uniform shifted 6x6x6 k-point grid is used for the 
continued fraction calculation. A 0.3 eV Lorentzian con- 
volution, gradually increasing to 1.5 eV after the edge is 
used to match the experimental width. 

III. RESULTS 

A. Experiment 

In figure [5] we show the measured PFY-XAS data at 
the Cu K-edge in C^-zCuC^C^ compared to previous 
PFY-XAS measurements [l[ on La2Cu04. In the case 
of Ca2~a;Cu02Cl2 the scattering geometries considered 
are: (i) k parallel to the ab plane and e parallel to the 
c-axis, labeled shortly e || c and (ii) k parallel to the c- 
axis and e parallel to the ab plane, labeled e || ab. In the 



case of La2Cu04_the geometries are slightly different, as 
reported in ref. [1|. For the e || ab geometry we also plot 
standard XAS data in Fig. [2] 

In the absence of a core-hole all d-states are filled ex- 
cept the 3d x 2_ y 2 orbital which is half-empty. However, 
the 3d a; 2_j ( 2 being planar, the quadrupolar part is for- 
bidden in both geometries. This is due to the fact that 
the quadrupolar matrix element of the absorption cross 
section includes a term (k- r)(e ■ r). In the e || c geometry 
e ■ r = 0, while in the e || ab geometry k • r = 0. So in 
both experimental settings the matrix element is zero. 

We first focus on the pre-edge region. At 8978 eV in 
e || c geometry an extremely small feature (labeled A in 
the following) is present, its intensity is 0.4 % of the edge 
jump. A similar feature is present at the same energy in 
the e || ab geometry (labeled A'), however the intensity 
in this direction is substantially larger, of the order of 3% 
of the edge jump. Interestingly this feature in the e || ab 
geometry merges with a second broader band centered 
at 8982 eV (labeled B'). Given the weakness of some of 
these excitations, particularly along the c-axis, we recon- 
sider our published PFY-XAS Cu-K-edge La2Cu04 data 
to see if similar features occur jjj . This is indeed the 
case, however these small features (in both geometries) 
were not discussed in our previous work llj. We have not 
found any XAS measurements d, 0, HE HH on cuprates 
displaying the A and A' features. We also performed 
XAS measurements on the most interesting e || ab direc- 
tion and found that the standard XAS spectrum shows 
a very broad tail extending down to the pre-edge region. 
None of the pre-edge features that are clearly resolved by 
PFY are visible here, as in the case of La2Cu04. Fur- 
thermore the B' at 8922 eV in the e || ab geometry is 
usually assumed to be completely quadrupolar. This is 
due to the lesser sensitivity of standard XAS data when 
compared with PFY measurements. 

In Ca2- x Cu02Cl2 the first feature is probably not 
more than 0.3 eV above the Fermi level and thus it is 
relevant for its low energy physics. The other pre-edge 
features at higher energies (labeled B, B' and C in fig. 
[3]) are common to those of other cuprates. 

The structure in the absorption edge is similar in the 
two systems, apart from a remarkable difference. In the 
case of Ca2-xCu02Cl2 with e || c geometry the edge 
(peak of the whiteline) at 8989 eV is accompanied by a 
distinct shoulder at 8992 eV. In Ca2~ x Cu02Cl2 for e || ab 
the edge seems substantially larger than the experimen- 
tal resolution, suggesting the presence of a second peak 
in this case too. In La2Cu04 there is no evidence of such 
a splitting. However the poorer experimental resolution 
we had in ref. [l] could hinder the detection of a second 
peak. Consequently higher resolution measurements are 
necessary to clarify if this splitting is a general feature of 
high T c cuprates. 
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FIG. 3: Comparison between calculated Cu K-edge 
XAS (solid lines) and experimental (dots) PFY-XAS of 
Ca2-a;Cu02Cl2 (x — 0.2) for both conditions of polarization 
e || c and e || ab respectively in the top and bottom panels. 
The core hole width of the calculation is 0.7 eV. The zero 
energy corresponds to 8979.2 eV 



B. Theoretical understanding 

1. pre-edge region 

To understand the origin of the features seen in the 
absorption spectra, we carried out extensive DFT cal- 
culations for Ca2_a;Cu02Cl2. In the presence of a core- 
hole the d x 2_ y 2 state around the absorbing atom is the 
only d state partially occupied. In the pre-edge region 
the situation is complicated by the fact that the sys- 
tem is metallic and it is hard to isolate the occupied 
states from the unoccupied ones, even as far as the dipo- 
lar part is concerned. As a consequence we have decided 
to shift the Fermi level 0.7 eV below that found in the 
supercell calculation (this is illustrated by the dashed 
line in Figs. I3I4I5|) . This shift improves the agreement 
with experiments in the low-energy region. However the 
drawback is that a small tail of the d xz + d yz + d z i pro- 
jected density of states becomes unoccupied generating 
a quadrupolar term (see Fig|5]). By direct calculation of 
the XAS quadrupolar cross section we have verified that 
this quadrupolar term is completely negligible with re- 
spect to the dipolar part emerging from the Fermi level 
downshift. Despite these difficulties and the strongly cor- 
related nature of the system which is a real challenge for 
density functional theory, we find good overall agreement 
with experiments. 

In the pre-edge region and e || c geometry, we repro- 



FIG. 4: Comparison between the calculated XAS cross sec- 
tion with c in the c direction and the density of states in the 
presence of a core-hole projected on selected orbitals. Note 
that in this figure the XAS has been plotted in log scale to 
enhance the low intensity low energy peaks. 



duce the low intensity peak A as it can be seen in fig. [3] 
The analysis based on the density of states projected on 
selected orbitals in the presence of a core-hole in the fi- 
nal state (fig. 2]) suggests that the peak is essentially due 
to an off-site transition to CI 2p z states and practically 
negligible quadrupolar component. This is analogous to 
what was found in NiO [2|, namely a very low energy 
dipolar feature due to direct off-site transition to nearby 
atoms. The situation is not as favorable in the e || ab ge- 
ometry where the broad band B' and the small feature A' 
are replaced in our calculation by a more intense broad 
band (figs. [3] and [5]). This disagreement is a reminder 
that DFT has limitations when it comes to the descrip- 
tion of the low energy physics of high T c compounds. 

On the other hand, the peaks C and D are correctly 
reproduced. The C peak is due to intersite hybridization 
between Cu 4p , O 2-p xy and Ca d states, while the D 
peak is due to transitions to on-site Cu 4p states. Sim- 
ilarly to the C peak, the B peak in e || c geometry is 
due to intersite Cu 4p, Ca d and CI p z states. Thus all 
the pre-edge peaks A B, and C are dominated by off-site 
transitions. 



2. Edge-to-far-edge region 

This energy region has been studied in detail in the 
past mainly due to the work of Tolentino et al. [f| and 
Kosugi et al. 0] suggesting that multi-determinant exci- 
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FIG. 5: Comparison between the calculated XAS cross section 
with e in the CuO plane direction and the density of states 
in the presence of a core-hole projected on selected orbitals. 



tations could occur. In particular, in ref. b] it was pro- 
posed that the edge peaks C and E' are due to transition 
to |3d 10 i) states while their supposed sattelites, E and 
F' are due to \3d?L). The original proposal was on Cu 
K-edge of La 2 Cu04. However in ref. § we showed that 
this system is not well described by DFT in this region 
due to the difficulty in treating La f states . Nevertheless 
we concluded in ref. Q that apart from a shift of the 
E peak to high energy, generated by an incorrect treat- 
ment of La f states in DFT, C,E,E', F' were most likely 
due to single particle excitations. The present compound 
is ideal to verify this conclusion since no rare-earths are 
present. Indeed we find that peaks C,E,E', F' are cor- 
rectly reproduced confirming their single particle nature. 

In the e || c direction a peak labeled D is found. This 
peak is absent from our single particle calculation. This 
suggests that multi-determinant charge transfer effects 
could indeed occur in cuprates. In order to test our con- 
jecture one should consider the convolution between hard 
X-ray Cu Is photoemission data and first principles cal- 
culations. This procedure should give semi-quantitative 
agreement with the Cu K-edge XAS [l9|-|3[. The diffi- 
culty in the present case is that there are no photoemis- 



sion data for this system so measurements and calcula- 
tions on other correlated compounds arc needed. 



IV. CONCLUSION 

In this work we have carried out Cu K-edge PFY-XAS 
measurements and state-of-the-art first principles calcu- 
lations in high T c cuprates Ca 2 _2;Cu02Cl2. We have 
shown that rich structure due to dipolar transitions oc- 
curs in the pre-edge region. By performing electronic 
structure calculations we have shown that pre-edge peaks 
are dominated by off-site transitions. Furthermore we 
have shown the occurrence of new excitations in the low 
energy region, at less than 0.5 eV from the Fermi level. In 
this energy region DFT calculations only partially cap- 
ture the relevant physics of the system since they fail 
in describing the in-plane excitations. Our work demon- 
strates the need of going beyond single site cluster models 
and, for what concerns correlation effects, beyond density 
functional theory, in order to explain Cu K-edge XAS of 
cuprates. 

Previous work invokes a strong influence of many 
body effects to explain the features observed in XAS K- 
edge spectra of high Tc materials. Notably the invoca- 
tion of a correlated ground state and the rearrangement 
of the energy states due to the core hole effect has lead 
earlier authors to conclude that some of the most promi- 
nent peaks in XAS spectra are not due to single parti- 
cle excitations. This interpretation is also inspired from 
core-level photoemission data but here we show that the 
situation is not analogous and that much of the XAS K- 
edge spectrum can in fact be explained invoking single 
particle excitations. Further, we find that many body 
effects due to charge transfer do occur but in a different 
way from earlier interpretations. 

The ramifications of this reinterpretation go beyond 
XAS spectroscopy. It will suffice to note that RIXS spec- 
tra use the resonances found in XAS to excite particular 
intermediate states. 
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